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Abstract. In a previous study it was proposed that the 
Galactic dark matter being detected by gravitational mi- 
crolensing experiments such as MACHO may reside in 
a population of dim halo globular clusters comprising 
mostly or entirely low-mass stars just above the hydrogen- 
burning limit. It was shown that, for the case of a stan- 
dard isothermal halo, the scenario is consistent not only 
with MACHO observations but also with cluster dynam- 
ical constraints and number-count limits imposed by 20 
Hubble Space Telescope (HST) fields. The present work 
extends the original study by considering the dependency 
of the results on halo model, and by increasing the sample 
of HST fields to 51 (including the Hubble Deep Field and 
Groth Strip fields). The model dependency of the results 
is tested using the same reference power-law halo models 
employed by the MACHO team. For the unclustered sce- 
nario HST counts imply a model-dependent halo fraction 
of at most 0.5 — 1.1% (95% confidence), well below the 
inferred MACHO fraction. For the cluster scenario all the 
halo models permit a range of cluster masses and radii to 
satisfy HST, MACHO and dynamical constraints. Whilst 
the strong HST limits on the unclustered scenario imply 
that at least 95% of halo stars must reside in clusters at 
present, this limit is weakened if the stars which have es- 
caped from clusters retain a degree of dumpiness in their 
distribution. 

Key words: stars: low-mass, brown dwarfs - globular 
clusters: general - Galaxy: halo - dark matter 



1. Introduction 

The abundance and nature of dark matter in the halo of 
our Galaxy is rapidly making the transition from theoret- 
ical hypothesis to observational science. This has been fa- 
cilitated by the deep surveys that are now achievable with 
instruments such as the Hubble Space Telescope (HST) 



and by several gravitational microlensing searches that are 
currently in progress. 

The 2nd-year results from the MACHO microlensing 
experiment (Alcock et al. 1997) towards the Large Mag- 
ellanic Cloud (LMC), a direction which is sensitive to 
lenses residing in the dark halo, indicate that a substan- 
tial fraction of the halo (20 — 100%) comprises objects 
with a typical mass in the range 0.1 — 1 Mq. 1 These re- 
sults appear to be broadly supported by the provisional 
findings from 4 years of MACHO observations (Axelrod 
1997) which have uncovered at least 14 LMC microlens- 
ing candidates. Similar mass scales have also been impli- 
cated by the EROS I microlensing experiment (Renault 
et al. 1997), though with a somewhat lower inferred halo 
fraction. These results are consistent with the lenses being 
in the form of low-mass hydrogen-burning stars or white- 
dwarf remnants. 

However, both of these candidates appear unattrac- 
tive when other observational and theoretical results are 
taken into consideration. The number density, age and 
mass function of white dwarfs in the Galaxy is strongly 
constrained by number counts of high-velocity dwarfs in 
the Solar neighbourhood, and by their helium production 
(Carr et al. 1984; Ryu et al. 1990; Adams & Laughlin 
1996; Chabrier et al. 1996; Graff et al. 1997). In particular, 
Chabrier et al. (1996) find that a halo fraction compatible 
with MACHO results requires that white dwarfs be older 
than 18 Gyr, though more recently Graff et al. (1997) have 
argued for a lower limit closer to 15.5 Gyr based upon rea- 
sonable white-dwarf model assumptions and a halo frac- 
tion of 30%. The situation for low-mass stars appears at 
least as pessimistic with recent HST results indicating that 
a smoothly distributed population of low-mass stars can 
contribute no more than a few percent to the halo dark 
matter density, regardless of stellar metallicity (Bahcall et 
al. 1994; Graff & Freese 1996; Flynn et al. 1996; Kerins 
1997). 
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1 This conclusion can be avoided if one instead attributes the 
observations to lenses residing in a very massive disc, though to 
explain the MACHO results one requires a local disc column 
density in excess of that typically inferred from kinematical 
observations (Kuijken & Gilmore 1989; Bahcall et al. 1992). 
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It has been suggested (Kerins 1997, hereafter Paper I) 
that if low-mass stars are clumped into globular-cluster 
configurations then HST limits can be considerably weak- 
ened, since this introduces large fluctuations in number 
counts and also may prevent a significant fraction of 
sources within the cores of clusters from being resolved. 
Motivation for the cluster scenario comes from the pre- 
dictions of some baryonic dark matter formation theories, 
which are discussed in Paper I. However, such clusters are 
required to have masses and radii consistent with exist- 
ing dynamical constraints on clusters and other massive 
objects residing in the halo. In Paper I it was shown that 
agreement between HST counts, dynamical limits and the 
central value for the halo fraction inferred by MACHO 
(40% for the halo model assumed) is possible if clusters 
have a mass around 4 x 10 4 Mq and radius of a few par- 
sees. However, HST, MACHO and dynamical limits are 
all dependent upon the unknown halo distribution func- 
tion, so these results are valid only for the spherically- 
symmetric, cored isothermal halo model adopted in Pa- 
per I. 

In this paper the model dependency of such conclu- 
sions is investigated using the same set of reference halo 
models employed in the MACHO collaboration's analy- 
sis of its results. One of these models is similar, though 
not identical, to the model investigated in Paper I, whilst 
the other models are constructed from the self-consistent 
set of 'power-law' halo models presented by Evans (1994). 
All models are normalised to be consistent with observa- 
tional constraints on the Galactic rotation curve and local 
column surface density. New data from the Hubble Deep 
Field (Flynn ct al. 1996) and Groth Strip (Gould et al. 
1997) are also incorporated, as well as two other new fields 
analysed by Gould et al., extending the analysis from 20 
HST fields in Paper I to 51 in this study. 

2. Halo models 

In Paper I constraints on the halo fraction in clustered and 
unclustered low-mass stars are derived assuming the stars 
have zero metallicity and that the halo density p varies 
with Galactocentric cylindrical coordinates (R, z) as 

where, in Paper I, the local density po = 0.01 Mq pc~ 3 , 
the Solar Galactocentric distance i? = 8 kpc, and the 
halo core radius R c = 5 kpc. 

The assumption of zero metallicity is maintained in 
the present analysis since one expects the halo to be per- 
haps the oldest of the Galactic components, and hence its 
constituents to have more or less primordial metallicity. 
The expected absolute magnitude in various photomet- 
ric bands for such stars between the hydrogen-burning 
limit mass (0.092 Mq) and 0.2 Mq has been calculated 



by Saumon et al. (1994) and their results are employed 
here as in Paper I. 

The model dependency of the conclusions in Paper I is 
assessed by re-calculating the constraints for a number of 
different, but plausible, halo models. For ease of compari- 
son the models selected are 5 of the reference halo models 
used by the MACHO collaboration (Alcock et al. 1996) 
in its analysis. (MACHO considers a total of 8 Galactic 
models, though only 5 of the halo models have distinct 
functional forms.) All halo models assume R = 8.5 kpc 
and R c — 5 kpc. The 5 models are denoted by MACHO 
as models A-D and S (for 'standard'), and this labelling 
is maintained here. 

The standard model S has the same functional form 
as the halo investigated in Paper I (i.e. it is described by 
Eq. 1) but uses the slightly larger IAU value for R above 
and assumes a lower local density po = 0.0079 Mq pc~ 3 . 
Models A-D are drawn from the self-consistent family of 
power-law models (Evans 1994), having density profiles 

_ vjRP R 2 J1 + 2g 2 ) + R 2 (l - /?g 2 ) + z 2 [2 - <T 2 (1 + (3)] 
9 AirGq {Rl + R 2 + z 2 q -2^/3+4)/2 

(2) 

where v a is the velocity normalisation, q describes the flat- 
tening of cquipotcntials, (3 determines the power-law slope 
of the density profile at large radii, and tt and G have their 
usual meanings. For a flat rotation curve at large radii 
(3 = 0, where as for a rising curve (5 < and for a falling 
one (3 > 0. 



Table 1. Parameter values for the 5 MACHO reference halo 
models A-D and S (Alcock et al. 1996). Rq = 8.5 kpc and 
R c = 5 kpc is assumed for all models. For models A-D the local 
density po is derived from the parameters in columns 2-4. The 
local rotation speed v is computed from the combined halo 
and disc mass within Rq. 



Model 


ii a /km s 1 


<7 


13 


Po/Mq pc 3 


vo/km s 


A 


200 


1 


0.0 


0.0115 


224 


B 


200 


1 


-0.2 


0.0145 


233 


C 


180 


1 


0.2 


0.0073 


203 


D 


200 


0.71 


0.0 


0.0190 


224 


S 








0.0079 


192 



The particular parameters for models A-D, along with 
those of model S are listed in Table 1. Model A is the 
closest analogy to model S within the power-law family 
of models, whilst model B has a rising rotation curve at 
large radii, model C a falling rotation curve, and model 
D a flattening equivalent to an E6 halo. When combined 
with the MACHO canonical Galactic disc (Alcock et al. 
1996), the models give values for the local Galactic rota- 
tion speed vq within 15% of the IAU standard value of 
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220 km s 1 and have rotation curves that are consistent 
with observations. 

3. HST observations and halo fraction constraints 

Gould et al. (1997) have calculated the disc luminosity 
function for M-dwarf stars using data from several HST 
WFC2 fields. These include 22 fields originally analysed 
by Gould et al. (1996), along with the Hubble Deep Field, 
28 overlapping fields comprising the Groth Strip, and 2 
other new fields: a total of 53 WFC2 fields. In Paper I, 
20 of the original 22 fields are analysed, the other 2 fields 
being omitted due to statistical problems introduced by 
their close proximity to some of the other fields (namely 
that clusters appearing in these fields could also appear in 
the other fields and thus be double counted) . In this study 
these 20 fields arc combined with the new fields analysed 
by Gould et al. (1997), making the total number of fields 
51. The nearest-neighbour separation between these fields 
is sufficiently large that double counting is not expected 
to be a problem for clusters of interest. (The overlapping 
Groth Strip fields are treated as a single large field for the 
purpose of this study.) 

The limiting and saturation /-band magnitudes for the 
fields are listed in Table 1 of Gould et al. (1997). The 
Groth Strip is treated as a single field with an angular 
coverage of 25.98 WFC2 fields (this accounts for overlaps) 
and magnitude limits corresponding to the modal values 
listed in Gould et al. (1997). As in Paper I, these limits 
are translated into star-mass dependent limiting distances 
by converting the line-of-sight extinction values listed in 
Burstein & Heiles (1984) to /-band reddenings and using 
the photometric predictions of Saumon et al. (1994) for 
zero-metallicity low-mass stars. The predictions for the V 
and / bands are well fit by the colour-magnitude relation 

/ = -11.45 (V - I) 2 + 40.7 (V - I) - 24.5 (3) 

for 1.27 < V-I < 1.57 (corresponding to 0.2 > m/M > 
0.092). 

The analyses for the unclustered and clustered scenar- 
ios proceed as in Paper I, except that the models listed in 
Table 1 of this paper now replace the model used there. 
The calculations for the cluster scenario, which are de- 
scribed in detail in Paper I, assume that the surface- 
brightness profiles of the clusters follow the King (1962) 
surface-brightness law and take into account cluster re- 
solvability, as well as line-of-sight overlap. 

Table 2 lists the results for the unclustered scenario. 
Within the 51 HST WFC2 fields analysed a total of 145 
candidate stars with 1.2 < V — I < 1.7 are found, imply- 
ing a 95% confidence level (CL) upper limit on the average 
number of 166 stars. This colour range spans the V — I 
colour predictions of Saumon et al. (1994) for stars with 
masses in the interval 0.092 — 0.2 Mq, where the lower 
value corresponds to the hydrogen-burning limit. Com- 
parison with the expected number tabulated in Table 2 



Table 2. Constraints on unclustered zero-metallicity low-mass 
stars in the Galactic halo arising from the detection of 145 
candidate stars within 51 HST WFC2 fields. The second and 
third columns give the expected number of detectable stars 
A^exp for a full halo (/h = 1) for stars with masses of 0.2 Mq 
and 0.092 Mq (the hydrogen-burning limit mass), respectively. 
The last two columns give the 95% confidence upper limit on 
the maximum halo fraction / max . 



^cxp /max 



Model 


0.2 Mq 


0.092 M 


0.2 Mq 


0.092 M 


A 


183 000 


24100 


9 x 10 -4 


0.007 


B 


248 000 


30 800 


6 x 10~ 4 


0.005 


C 


109 000 


15 000 


0.0015 


0.011 


D 


162 000 


34300 


0.0010 


0.005 


S 


141 000 


17100 


0.0012 


0.010 



clearly shows that, for all models, even the lowest mass 
unclustered stars fall well short of providing the halo dark 
matter density inferred by MACHO. The upper limit on 
their fractional contribution / max is shown for 0.2-Mq and 
0.092-Mq stars. For the lowest mass stars / max ranges 
from 0.5% for models B and D to 1.1% for the lighter halo 
model C. 

One interesting feature of Table 2 is that for the flat- 
tened halo model D the expected number counts are en- 
hanced for 0.092-Mq stars relative to the predictions for 
the spherically symmetric models, producing the high- 
est predicted number-count for these stars. This contrasts 
with the results for the brighter 0.2-Mq stars, with the 
heavy halo model B producing the highest number-count 
prediction. The enhancement for 0.092-Mq stars in model 
D arises because the flattening preferentially increases the 
stellar surface density near the Galactic plane, and this 
is reflected in the counts of 0.092-Mq stars which can be 
at most only a few kpc from the plane if they are to be 
detected. 

The constraints on the halo fraction /h for the clus- 
tered scenario as a function of cluster mass M and radius 
R are shown in Fig. 1 for the 5 models (A-D, S) assuming 
all stars reside in clusters and have the hydrogen-burning 
limit mass of 0.092 Mq. Each plot is characterised by a 
lower plateau to the left, an upper plateau to the right and 
a curved rising surface joining the two. This curved sur- 
face between the two flat regions represents the 95% CL 
upper limit halo fraction in clusters inferred from the pres- 
ence of only 145 candidate stars within the 51 HST WFC2 
fields. The constraints are actually calculated on the basis 
of no stars being present within these fields, since for clus- 
ters there is little difference in the constraints assuming 
no stars are found or assuming a few hundred stars are 
found. The reason for this, as discussed in Paper I, is that 
the clusters considered here contain between 1000 and 10 7 
members each, so the presence of just one cluster within 
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Fig. 1. Comparison of constraints on the halo fraction /h from HST limits, MACHO observations and dynamical constraints 
for the 5 halo reference models (A-D, S), assuming halo stars have a mass of 0.092 Mq and all reside in clusters with mass 
M and radius R. The lower plateau to the left of each plot corresponds to the 95% CL upper limit /max for the unclustered 
scenario inferred from HST counts (see Table 2). The upper plateau on the right corresponds to the 95% CL lower limit halo 
fraction /m,1ow inferred by MACHO 1st- and 2nd-year observations (Alcock et al. 1997), with the central value for the MACHO 
halo fraction /m indicated by the skirting surrounding the plots (see also Table 3). The curved surface joining the lower and 
upper flat regions corresponds to the 95% CL upper limit on the halo fraction in clusters from HST counts. Also projected onto 
the plane /h = /m,1ow are the cluster dynamical constraints (dashed lines) for the local Solar neighbourhood. The intersection 
between these constraints and the MACHO lower-limit plateau indicates cluster parameters compatible with HST, MACHO 
and dynamical constraints. 
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any of the HST fields would typically result in thousands 
if not millions of candidates being detected. 

The lower plateau shows the 95% CL upper limit halo 
fraction for the unclustered scenario (corresponding to the 
/max values listed in Table 2). Clusters with masses and 
radii within this region have internal densities which are 
lower than that of the halo background average and are 
thus unphysical, since they represent local under-densities 
rather than over-densities. Clearly constraints on clusters 
cannot be stronger than constraints on a smooth stellar 
distribution. The intersection of the lower plateau with 
the curved rising surface therefore denotes the boundary 
between unphysical and physical cluster parameters. 



Fig. 1. continued. 



The upper plateau to the right represents the 95% CL 
lower limit on the halo fraction /m,1ow inferred from MA- 
CHO 1st- and 2nd-year microlensing results (Alcock et al. 
1997). It is calculated by taking the 95% CL lower limit on 
the measured microlensing optical depth for all 8 MACHO 
events (t > 1.47 x 10~ 7 ), subtracting the optical depth 
contribution expected from non-halo components [corre- 
sponding to T non _haio — 5 x 1CP 8 (Alcock ct al. 1996)], 
and normalising to the optical depth prediction r cxp for a 
full halo (/h = 1) for each model. The top of the skirting 
surrounding each plot is normalised to the central MA- 
CHO value for the halo fraction /m for comparison, and 
is calculated in a similar manner to the lower limit (/m 
and /m,1ow, together with T cxp , are tabulated in Table 3 
for each model). Since this plateau lies below the extrap- 
olation of the HST cluster-fraction constraint [which rises 



asymptotically over this region - see Fig. 2 of Paper I], it 
is consistent with both MACHO and HST limits. 

The dashed lines in the plots of Fig. 1 represent the dy- 
namical constraints derived for the local Solar neighbour- 
hood. In fact some of the HST fields are somewhat closer 
in to the Galactic centre, where the dynamical constraints 
are stronger, but most are further away so the limits shown 
are stronger than applicable for most of the HST fields. 
The functional form for the constraints are detailed in 
Paper I and are dependent upon Galactic as well as clus- 
ter parameters [consult Lacey & Ostriker (1985); Carr & 
Lacey (1987); Moore (1993); Moore & Silk (1995); Carr & 
Sakellariadou (1997) for derivations, and see Carr (1994) 
for a detailed review of dynamical constraints] . Their vari- 
ation from plot to plot is due to model variations in the lo- 
cal density and rotation speed (see Table 1). The dynam- 
ical constraints are projected onto the plane /h = /m,1ow 
for direct comparison with the MACHO lower limits. The 
intersection of the MACHO lower-limit plateau with the 
dynamical limits therefore represents cluster parameters 
compatible with MACHO, dynamical limits, and the con- 
straints from the 51 HST fields. 

For each model it is evident that the region compati- 
ble with all limits spans a significant range of masses and 
radii. For models C and D the maximum permitted clus- 
ter mass is around 5 x 10 4 Mg, whilst for models A, B 
and S one can have cluster masses in excess of 10 6 Mg. 
Interestingly, whilst in the unclustered scenario the heavy 
halo model B is the most strongly constrained in terms 
of allowed halo fraction / max , it nonetheless allows a rela- 
tively wide range of viable cluster masses in the clustered 
scenario. Conversely, the permitted cluster mass range for 
the light halo model C is more restricted. 

This apparent paradox is due to the fact that the HST, 
dynamical and microlensing observations limit the halo 
density normalisation at different positions in the halo, so 
their intersection is sensitive to the halo density profile. In 
particular, the HST and dynamical limits essentially apply 
to the local Solar neighbourhood position (R = 8 kpc) 
for clusters comprising relatively dim hydrogen-burning 
limit stars, where as the microlensing observations to- 
wards the LMC constrain the density of lenses at some- 
what larger distances (primarily between 10 and 30 kpc 
from the Galactic centre, where the product of lens num- 
ber density and lensing cross-section is largest). Hence, 
for a given microlensing constraint on the mass density of 
lenses at 10 to 30 kpc, the local dynamical and number- 
count constraints are weaker for haloes with rising rota- 
tion curves (such as model B) than for models with falling 
rotation curves (such as model C). 

The relatively large range in allowed cluster masses 
and radii for model S is in apparent contrast to the re- 
sults of Paper I, in which the surviving parameter space 
is shown to be much smaller for the very similar model 
adopted there. There are two reasons for this apparent 
discrepancy: (1) in Fig. 1 of this paper it is assumed that 
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the clusters comprise hydrogen-burning limit stars, where 
as in Fig. 3 of Paper I the constraints are shown for the 
brighter 0.2-Mq stars; (2) in this study consistency is be- 
ing demanded only with the lower limit MACHO halo 
fraction /m,1owi rather than with the central value /m as 
in Paper I. This latter difference is particularly impor- 
tant because it enlarges both the sizes of the dynamically- 
permitted region and the MACHO plateau, and hence en- 
larges their intersection. Since these differences serve to 
maximise the size of the surviving region, the constraints 
shown in this paper should be taken as firm limits on al- 
lowed cluster parameters. 

4. Constraints on cluster membership 

Figure f assumes that all stars reside in clusters at the 
present day, an unrealistic assumption since one expects 
some fraction of the clusters to have evaporated away over 
time. As in Paper I one can place limits on the fraction of 
stars f c which must remain in clusters by using the strong 
limits /max on the unclustered scenario (listed in Table 2). 
Assuming the lower limit on the cluster halo fraction to be 
given by the lower limit inferred by MACHO, /m.Iow, the 
present-day halo fraction in stars which have evaporated 
away from clusters is fh,* > (f — / c )/m,1ow Since HST 
observations demand fh,* < /max one has 

fc > 1 - (/max//M,low)- (4) 

The resulting values for f c for 0.2-Mq and 0.092-Mq stars 
are given in Table 3. 



Table 3. Microlensing halo fractions and minimum cluster- 
ing fractions for the reference halo models. Column 2 gives 
the expected optical depth for a full halo as calculated by Al- 
cock et al. (1996). Column 3 gives the central value for the 
halo fraction using the lst+2nd year optical depth estimate of 
2.94 x 10 -7 measured by Alcock et al. (1997), and subtracting 
from it an optical depth of 5 x 10~ 8 expected from non-halo 
populations. The 4th column gives the 95% CL lower limit on 
the halo fraction using the lower limit for the measured optical 
depth of 1.47 x 1CP 7 . The last two columns give the lower limit 
on the present-day clustering fraction using column 4, Eq. 4 
and Table 2. 



fc 



Model 


T Cxp /i(n 7 


/m 


/m,1ow 


0.2 M 


0.092 M 


A 


5.6 


0.43 


0.17 


0.995 


0.96 


B 


8.1 


0.30 


0.12 


0.995 


0.96 


C 


3.0 


0.81 


0.32 


0.995 


0.97 


D 


6.0 


0.41 


0.16 


0.994 


0.97 


S 


4.7 


0.52 


0.21 


0.994 


0.95 



From Table 3 it is clear that all models require a very 
high fraction of all stars to reside in clusters at present. 



Even for hydrogen-burning limit stars the required clus- 
tering fraction must be at least 95% at present. Capriotti 
& Hawley (1996) have undertaken a detailed analysis of 
cluster mass loss within an isothermal halo potential for 
a range of cluster masses, density profiles and Galacto- 
centric distances. Their analysis takes account of evapora- 
tion, disruption and tidal processes. They find that clus- 
ters with masses between 10 5 — 10 7 Mq generally survive 
largely intact to the present day but that less massive clus- 
ters survive only if they have high central density concen- 
trations, nearly circular orbits and reside at large distances 
from the Galactic centre. At the Solar position Capriotti 
& Hawley find that clusters with a half-mass to tidal radii 
ratio of 0.3 (comparable to the value for the clusters anal- 
ysed here and in Paper 1) survive more than 95% intact 
only if they have masses exceeding 10 6 Mq. 

However, there are a number of reasons why these lim- 
its may be stronger than applicable to the low-mass star 
cluster scenario. Firstly, Capriotti & Hawley assume that 
the clusters comprise m = 0.8 Mq stars (i.e. between 4 
and 9 times more massive than the stars considered in the 
present study). The evaporation timescale scales approx- 
imately as mT 1 for a fixed cluster mass, so for clusters 
comprising lower mass stars the evaporation timescale is 
correspondingly longer. Secondly, the local halo density 
assumed by Capriotti & Hawley of 0.0138 Mq pc~ 3 at a 
Galactocentric distance of 8.5 kpc is on the higher end 
of the values for the halo models analysed in this paper, 
and is considerably larger than the allowed MACHO lower 
limit, /m,1owP0j ° n the local density in lenses (by a fac- 
tor of between 5 and 10 after normalising to a distance 
i?o = 8 kpc). Hence disruption due to close encounters 
with other clusters is substantially less in the halo models 
investigated here than for the model analysed by Capriotti 
& Hawley. 

Lastly, a study by Oh & Lin (1992) has shown that 
the cluster escape rates may be substantially smaller than 
commonly assumed due to angular momentum transfer 
arising from the action of the Galactic tidal torque on 
cluster stars with highly eccentric orbits (which in the 
absence of the torque would constitute the bulk of the 
escapees). The rates calculated by Oh & Lin for isotropic 
cluster models are broadly consistent with the values used 
by Capriotti & Hawley (1996) and other authors, but for 
the case of anisotropic stellar orbits the escape rates can 
be 1-2 orders of magnitude smaller, again implying cor- 
respondingly longer evaporation timescales. ft therefore 
appears that, under certain conditions, one may be able 
to reconcile the high cluster fraction requirements derived 
in the present study with the findings of cluster dynamical 
studies, at least for clusters comprising stars close to the 
hydrogen-burning limit. 

In any case, the validity of the figures in Table 3 de- 
pend upon just how smoothly distributed are the stars 
which have evaporated from clusters. If they still have 
not completely homogenised today, instead maintaining a 



E.J. Kerins: Low-mass stars and star clusters in the dark Galactic halo 



7 



somewhat lumpy distribution (reflecting their cluster ori- 
gin), then the limits on f c are too strong. 

For example, a cluster with a mass 3 x 10 4 Mq and ra- 
dius 3 pc represents an over-density of about 3 x 10 4 over 
the background average at the Solar neighbourhood [i.e 
Sp/p = (p — p)/p = 3 x 10 4 ]. However, an under-density in 
the unclustered (or more precisely 'post-clustered') stellar 
population of just a factor 10 (Sp/p = —0.9) over vol- 
umes larger than 3 x 10 5 pc 3 , which is roughly the volume 
probed by 50 HST fields for hydrogen-burning limit stars 
(and is of order 10 times smaller than the halo volume per 
cluster) , is all that is required to weaken the constraints on 
/max by a factor 10. This would result in a much more com- 
fortable lower limit on f c of just 0.5 for 0.092-Mq stars, 
rather than 0.95. If the under-density is a factor 5 lower 
than the background (Sp/p = —0.8) one requires / c > 0.75 
for the lowest mass stars and for an under-density factor 
of 2 (Sp/p = —0.5) f c must exceed 0.9. 

In order to rule out the cluster scenario definitively 
(say with 95% confidence) one needs a survey that is both 
sufficiently wide and deep that it might be expected to 
contain at least 3 clusters on average, regardless of their 
mass and radius (though their mass and radius must be 
dynamically permitted). From Fig. 1 it appears that the 
most difficult dynamically-allowed clusters for HST to ex- 
clude arc those with a mass of around 3 x 10 4 M m . If the 
halo fraction in low-mass stars is around 40%, typical of 
the preferred value for the MACHO results, then the lo- 
cal number density of such clusters is around 130 kpc~ 3 
(adopting a local halo density of 0.01 Mq pc~ 3 ; in reality 
of course the average density within the fields is depen- 
dent upon the halo model and the field locations). If the 
clusters comprise hydrogen-burning limit zero-metallicity 
stars (V — I = 1.57) then a HST-type survey will be sensi- 
tive to them out to about 3.6 kpc in the I band [using the 
colour-magnitude relation of Eq. (3), and assuming a lim- 
iting /-band sensitivity of 24 mag] , so in order to expect 
to detect at least 3 such clusters, the survey must cover a 
solid angle of at least 4.5 deg 2 , or the equivalent of 3 700 
HST fields! Therefore, only if HST fails to detect any clus- 
ters from 3 700 fields would dynamically-allowed clusters 
be ruled out with 95% confidence from explaining all of the 
observed microlensing events. For comparison, an all-sky 
X-band survey over Galactic latitudes \b\ > 10° requires a 
limiting magnitude of about 17.5 in order to produce sim- 
ilar constraints. This should be compared to the expected 
if-band limit of about 14 for the ground-based DENIS 
and 2MASS surveys. 

An easier alternative is to instead obtain several fields 
as close to the Galactic centre as is feasible, where the 
dynamical constraints are much stronger than for the So- 
lar neighbourhood position. For low-mass stars, this ne- 
cessitates a telescope such as HST with the capability 
for obtaining very deep fields, since shallow surveys with 
wide angular coverage essentially only probe the local So- 
lar neighbourhood. 



5. Conclusion 

Kerins (1997), referred to as Paper I, has suggested that 
low mass stars could provide the substantial dark mat- 
ter fraction indicated by the combined 1st- and 2nd-year 
MACHO gravitational microlensing results. Whilst obser- 
vations from Hubble Space Telescope (HST) and other 
instruments have been interpreted as excluding such stars 
from having a significant halo density, Paper I shows that 
their density could in fact be substantial if they are 
grouped into globular-cluster configurations. The motiva- 
tion for such clusters comes from the baryonic dark matter 
formation scenarios which are discussed in Paper I. Pa- 
per I calculates the constraints on such clusters (assum- 
ing they comprise low-mass stars of primordial metallicity) 
which arise from MACHO microlensing results, dynami- 
cal constraints on massive halo objects, and observations 
from 20 HST fields obtained by Gould et al. (1996). How- 
ever, the results of Paper I apply only to the spherically- 
symmetric cored isothermal halo model investigated there. 

In the present study, the number of HST fields utilised 
has been increased to 51, and now incorporates the Hub- 
ble Deep Field and Groth Strip fields (Gould et al. 1997). 
The model dependency of the results in Paper I has been 
tested by adopting 5 of the reference halo models employed 
in the MACHO collaboration's analysis of its microlens- 
ing results. One of the models is similar to the halo in- 
vestigated in Paper I whilst the other 4 are drawn from a 
self-consistent family of power-law halo models and com- 
prise spherically-symmetric haloes with a rising rotation 
curve, a falling rotation curve and a flat curve, as well as 
a flattened (E6) halo model. 

The 51 HST fields contain just 145 candidates with 
V — I colours between 1.2 and 1.7 (spanning the colour 
range predicted for zero-metallicity stars with masses be- 
tween the hydrogen-burning limit and 0.2 Mq) against 
the tens or hundreds of thousands predicted for the halo 
models. From this one concludes that the halo fraction 
in unclustered low-mass stars is at most 0.5 — 1.1% with 
95% confidence, depending on the halo model, and in all 
cases falls well short of providing even the lower-limit halo 
fraction inferred by MACHO. 

However, in the cluster scenario there exists a wide 
range of cluster masses and radii which can allow a halo 
fraction consistent with the lower limit derived from MA- 
CHO microlensing results whilst remaining compatible 
with dynamical limits and HST observations. Consistency 
with the preferred microlensing halo fraction, rather than 
the lower limit, requires fine tuning of the cluster param- 
eters (as found in Paper I), but is possible for all models 
investigated. 

The one potentially serious problem for the cluster sce- 
nario is that the strong constraints on unclustered stars 
imply that an overwhelming fraction of all stars, at least 
95%, must still reside in clusters at the present day. This 
is higher than expected from generic cluster evaporation 



8 



E.J. Kerins: Low-mass stars and star clusters in the dark Galactic halo 



considerations for much of the permitted cluster mass 
range, though it may still be consistent with clusters com- 
prising stars with anisotropic orbits. In any case, these 
limits assume that stars which have already evaporated 
from clusters now form a perfectly smooth distribution 
which traces the halo density profile. If instead these stars 
still have a lumpy distribution, reflecting the fact that 
they previously resided in clusters, then the cluster frac- 
tion limits are too strong. 

Probably the only way to definitively exclude or con- 
firm the cluster scenario is to obtain several deep fields 
as close to the Galactic centre as is practical, where the 
strong dynamical constraints severely restrict the range of 
feasible cluster parameters. 
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